Modified fluctuation-dissipation and Einstein relation at non-equilibrium steady states 



o 

(N 
S3 

(N 
(N 



^3 

o 



c3 



S3 
O 

o 



in 

(N 

in 
o 

(N 



X 



Debasish Chaudhuri^B and Abhishek Chaudhuri 2 ' 3 '0 

1 FOM Institute for Atomic and Molecular Physics, 
Science Park 104, 1098XG Amsterdam, The Netherlands^ 
2 Department of Biomedical Science, University of Sheffield, 
Western Bank, Sheffield S10 2TN, United Kingdom 
3 The Rudolf Peierls Centre for Theoretical Physics, 
University of Oxford, 1 Keble Road, Oxford 0X1 3NP, United Kingdom 
(Dated: January 24, 2012) 

Starting from the pioneering work of G. S. Agarwal [Zeitschrift fur Physik 252, 25 (1972)], we 
present a unified derivation of a number of modified fluctuation-dissipation relations (MFDR) that 
relate response to small perturbations around non-equilibrium steady states to steady-state corre- 
lations. Using this formalism we show the equivalence of velocity forms of MFDR derived using 
continuum Langevin and discrete master equation dynamics. The resulting additive correction to 
the Einstein relation is exemplified using a flashing ratchet model of molecular motors. 
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I. INTRODUCTION 

Derived within linear response theory, the fluctuation 
dissipation theorem (FDT) predicts how the response 
function of a thermodynamic observable is related to cor- 
relation of thermal fluctuations at equilibrium. Let us as- 
sume that an equilibrium system described by a Hamil- 
tonian H is perturbed at time t = t\ by an external 
force h(t). The FDT predicts a response at a later time 
ts > h Q 



R7(h h) = = pJL {A{t2)hdhH{tl)]h=o)eq 

(i) 

where the correlation is calculated at equilibrium cor- 
responding to temperature T with ft = 1/T. The dif- 
ferential operator dh in the above relation denotes the 
scalar derivative evaluated at time ti. Thus —dhH is the 
displacement conjugate to h with respect to the Hamilto- 
nian. Throughout this paper we use Boltzman constant 
ks = 1, unless otherwise stated. Using the Onsager re- 
gression hypothesis the FDT can be interpreted as fol- 
lows - the decay of a fluctuation is independent of how it 
has been created, under the influence of a small applied 
force or spontaneously by thermal noise. The FDT is vi- 
olated away from equilibrium regime and this violation 
has been studied in context of glassy systems, granular 
matter, sheared fluid, stochastic processes, and biological 
systems [2-|{|. 

In a pioneering study back in 1972 [l(|, G. S. Agar- 
wal obtained a modified fluctuation-dissipation relation 
(MFDR) that related response functions around non- 
equilibrium steady states (NESS) to correlations evalu- 
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ated at steady state. For a system evolving with a statis- 
tical dynamics characterized by the Fokker-Planck (FP) 
equation d t p = C$p, Agarwal showed that a perturbation 
in the operator Cq — > Co + h(t)Ci leads to a response that 
can be expressed in terms of a correlation function eval- 
uated at the unperturbed steady state [1(3, EH , 



R A (h-h) 



5(A(h)) 
Sh{h) 



(2) 



where the Agarwal term M = [Cip s ]/p s with p s denot- 
ing the steady-state probability distribution. Through- 
out this paper by (. . .) we denote a steady-state average. 

Over the last decade a formalism of stochastic thermo- 
dynamics has been developed that allows description of 
energy and entropy along fluctutating trajectories [l2| — 
fl~H . Various fluctuation theorems involving the distribu- 
tion of entropy (TB-dH , and work theorems [22T - |24| were 
discovered. Recently, using an integral fluctuation theo- 
rem, a number of these relations were derived in a uni- 
fied manner [lil , |25| . Important experimental tests in- 
clude colloidal particles manipulated by laser tra ps 1261- 
, biomolecules pulled by AFM or laser tweezer 29L l30j] 
and autonomous motion of motor proteins [3l| . Stochas- 
tic thermodynamics has also been used to derive several 
versions of MFDR around NESS @, 0, H Some 
of these predictions were experimentally verified (27l . [37| . 

In this paper, we present a unified derivation of a num- 
ber of MFDRs based entirely on the Agarwal formal- 
ism [l(|. Thus the MFDRs we obtain are intrinsically 
equivalent to each other. We show that the Agarwal 
term M can be expressed as a velocity excess from a 
local mean velocity using both the continuum Langevin 
and discrete master equation dynamics. This interpreta- 
tion leads us to a modified Einstein relation that has the 
same additive correction term for the two cases. Finally 
we apply this framework to a flashing ratchet model of 
molecular motors I3&f40l to calculate the MFDR and the 



additive correction in Einstein relation, which shows a 



non-monotonic variation with the asymmetry parameter 
of the ratchet. 

The structure of this paper is as follows. In Sec. HT1 we 
review the derivation of the Agarwal form of MFDR, that 
we use throughout this paper to calculate other versions 
of MFDR expressed in physically observable form. Using 
this result, in Sec. IIIII we present a simple and straight- 
forward derivation of the MFDR in terms of stochastic 
entropy production, keeping in mind that this relation 
was used earlier to derive velocity-MFDR for a master 
equation dynamics @ . Then, directly using the Agarwal 
form, we derive the velocity-MFDR for a system evolv- 
ing with continuum Langcvin dynamics in Sec. IIV1 and a 
discrete master equation in Sec. [V] The velocity-MFDR 
is used in Sec. I VII to derive a modified Einstein relation 
at NESS. In SecHU we study the velocity-MFDR, and 
the violation of the Einstein relation in a flashing ratchet 
model of molecular motors. Finally in Sec. I Villi we sum- 
marise our main results and conclude. 



probability. The time evolution dtp = Cop can be 
solved to obtain the transition probability at steady state 
iu(?, t|<r', 0) = exp(£ot)5(<; — Thus the two-time cor- 
relation at steady state takes the form (A(t)B(0)) = 
J dsA(s) exp(£ t)B((;)p s (<;). Therefore we can write 
Eq. [6] as 

R A (t2-t 1 ) = {A(t 2 )M(t 1 )). (7) 

This is the Agarwal form of MFDR [lo|. The derivation 
presented here used a continuum notation of the phase 
space variable However, the result is general, and can 
be derived similarly for a system that evolves through 
transitions between discrete states (see Eq. l2Tj) . 

The Agarwal term in its operator form M(s) = 
[Cip s ]/p s requires detailed knowledge of the probability 
distribution at steady state. In the rest of this paper 
we focus on expressing this term in physically observable 
form. 



II. THE AGARWAL FORM OF MFDR 



III. MFDR IN TERMS OF STOCHASTIC 
ENTROPY 



The probability distribution p(s, t) of finding a system 
at state <r at time t evolves with time as 



d t p(<;,t) = C(s,h)p(<;,t) 



(3) 



where £ is a general time evolution operator that depends 
on external force h(t). For weak h, Taylor expanding the 
operator we get 

C{s,h) = £o{s) + h{t)Cxk) (4) 
where C\ = [dh.C]h=o- The solution to Eq. [3]is 

p(s,t)=p 3 + [ dTe Co ^h(T)C lPs ^) (5) 



where p s denotes the steady-state distribution that obeys 
CoPs = 0. Then the response of any observable (A(t)) = 
J dsA((;)p(q,i) to a force h(t) is 



$p(<;,t 2 ) 

Sh{h) J Sh(h) 
= f d^A^e^-^dp^) 



J dsA^e^-^M&p^) 



(6) 



where in the last step we used the Agarwal term M(cr) = 
[Cip s ]/p s . By definition, the two-time correlation func- 
tion is (A(t)B(0)) = jdsjdq>A{<;)B{<;')p 2 {<;,tW,Q), 
where p 2 (<;,t\ 0) is the joint probability distribution 
of finding the system at state at time and at 
state ? at time t. One can express p 2 (s,t;<;',0) — 
«;(?, t|?',0)p(?',0) where w(<;, t\q' , 0) is the transition 



The definition of non-equilibrium Gibb's entropy S = 

— J dsp(s,t) lnp(<r, f) = (s(i)) has recently been used 
to get a definition of the stochastic entropy s(t) = 

— lnp(<r,t) [25j]. For a master equation based discrete 
dynamics between states denoted by n(t), the stochas- 
tic entropy can be written as s(t) = — ]np n n\. Using 
this definition we obtain a simple interpretation of the 
Agarwal term in terms of stochastic entropy 



M 



— ClPs 
Ps 

dhdtP 



d h C(h)p 



h=0 



p 



-d t [d h s\ 



h=0- 



(8) 



h=0 



In deriving the above relation we assumed that C{h) is 
linear in h. We also used the fact that the steady state 
distribution p s = p\h=o- Thus M is expressed as time- 
evolution of a variable conjugate to the external force h 
with respect to the stochastic system-entropy s. In this 
sense, s in NESS plays the role similar to the Hamiltonian 
in equilibrium FDT. We can now write the MFDR at 
NESS as 

R A (t 2 -h) = ^-(A(t 2 )[-d h s(h)} h = ). (9) 
at i 

Ref. [3, EH found this relation by considering a perturba- 
tion that takes the system to a final steady state. Note 
that our simple and straightforward derivation does not 
require such an assumption, and thus the result is more 
general. 



A. Equilibrium FDT 

The FDT at equilibrium can easily be derived from 
Eq. [5] If, even in the presence of external perturbation 



2 



the system remains at equilibrium, one can write down 
the probability distributions as p = cxp[—f3(H — F)} 
where F is the free energy. This distribution leads to 
the relation [d h p\ h =a = /3[(d h F - d h H)p) h=0 . Note 
that the equilibrium displacement evalutaed at h = 
is [di l F]i l= o = 0. Thus we get the identity [dhs]h=o = 
— [(9hP)/p]h=o = P[dhH]h=o, which leads to the equilib- 
rium FDT Eq. [TJ 



IV. VELOCITY MFDR USING LANGEVIN 
EQUATION 

Let us consider a Langcvin system where the dynamics 
of a particle evolves by 



v = (if + rj 



(10) 



where v = x is the particle velocity, fi is the mobil- 
ity, and / denotes total force imparted on the particle. 
The total force f(x,t) consists of a force due to inter- 
action F{x) and an external time dependent force h(t): 
f(x,t) = F(x) + h(t). The last term 77 denotes a ther- 
mal noise that obeys (77) = and (77(4)77(0)) = 2DS(t) 
with D = fiT, the equilibrium Einstein relation. The 
corresponding FP equation is 



d t p(x, t) 
with, j{x,t) 



-d x j{x,t) (11) 

(fif(x,t)-Dd x )p(x,t). 



The velocity form of MFDR for a Langevin system was 
originally derived in Ref . [lij . Here we briefly outline the 
derivation starting from the Agarwal form. Eq. [11] can 

be expressed as, 

d t p{x, t) = (C Q + h(t)d)p(x, t), 

where, Co = —d x (/iF) + Dd x , and C\ = —fid x . Thus 
the Agarwal term M = —fi(d x p s )/p s , and TM = 
—D(d x p s )/p s . The definition of the steady state current 
j s leads to the relation Dd x p s = /iF(x)p s (x) —j s . Defin- 
ing a local mean velocity at steady state v s (x) — j s /p s (x) 
we can then rewrite TM = —D(d x p s )/p a = v s (x) — 
fiF(x). In this relation, using the Langevin equation at 
initial steady state (h = 0), we get TM = v s — v + 77. 
Thus, the response function 

TR A (t 2 - ti) = (A(t 2 )[v(ti) - v(h) + ij(ti)]). (12) 

Note that in the Langevin equation /j,h(t) and 77(7;) have 
the same status, and A(x, t) can be regarded as a func- 
tional of noise history. Then it can be shown that [34| , 

TR A (t 2 h) = D 5 -^^- = ±{A(hM tl )). (13) 

Thus we can write Eq. [12] as 

R A (t 2 -h)=p (A(t 2 ) [v(h) - (14) 



This is the velocity form of MFDR, which for velocity- 
response gives 



Rv{h -h) = p (v(h)[v{h) - z/(*i)]>. 



(15) 



Note that the steady state average of v is the same as 
the mean velocity: 



L/2 



-L/2 



dxp s (x)v s (x) = /j,(F) - D\p s ]^ /2 = (vg) 



(16) 



The boundary term [p s ]^ 2 y 2 = either by a periodic 
boundary condition (34|, or by taking the boundaries to 
infinity where the probabilities vanish. If the system is 
at equilibrium v = 0, and we get back the well-known 
equilibrium response, 



R e v q (t 2 ~h) = HHhHh)) eq . 



(17) 



Therefore the non-equlibrium MFDR Eq. [15] can be 
viewed as the equilibrium FDT with an additive correc- 
tion -P(v(t 2 )v(ti)). 

It is interesting to note that using Eq. 1131 we can 
arrive at a non-equilibrium MFDR first obtained in Ref. 
[32I ] for continuous Langevin dynamics and subsequently 
shown to be true for discrete spin variables (as well as, for 
conserved and non-conserved order parameter dynamics) 
in Ref. [33| . Defining the position correlation function 
C x (t 2 ,ti) = (x(t2)x(t\)) and the corresponding response 
function 2TR x (t 2 ,t 1 ) = (aj(i 2 )77(*i)) (using Eq. we 
get the modified MFDR 

(d tl -dt^Cxfah) = 2TR x (t 2 ,t 1 )+A(t 2 ,h) (18) 

where A(t 2 ,ti) = (fif(ti)x(t 2 ) - nf(t 2 )x(t{)) is the 
so-called asymmetry which vanishes in the presence of 
time reversal symmetry. Note that causality demands 
that the response of the system at time t 2 to a per- 
turbation at time ti, R x (t 2 ,ti), is nonzero only when 
t 2 > ti. Incorporating time translation invariace and 
time reversal symmetry restores the equilibrium FDT, 
TR x (t 2 ,ti) = du C x (t 2 ,U )■ Also note that, the choice of 
the observable V in Ref. [3y| as an ID coordinate x, re- 
duces the second term on the r.h.s of Eq. 13 in Ref. [3(| 
to ((L-L*)V(s)Q(t)) = (2(j/p)VxQ(t)) = (mQ). Now 
setting Q = v (velocity), leads Eq. 13 in Ref. [36| to Eq. 
[14] in our manuscript, the velocity form of MFDR. 



V. VELOCITY MFDR USING MASTER 
EQUATION 

We now focus on a master equation system where the 
time-evolution occurs via transitions between discrete 
states. Following Ref. [34j . we first derive the discrete 
form of the Agarwal term M. Our main contribution in 
this section is to express M as an excess velocity, and 
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thus arrive at a velocity form of MFDR, similar to the 
Langevin system. 

We begin by considering a set of discrete states {n} 
and write down the corresponding master equation for 
the probability p m (t) of finding the system in a state m 
at time t: 

d t Pm(t) = y^\w nm p n (t) - w mn p m (t)} 

n 

= Y, C rnnPn(t) (19) 
n 

where, w mn represents transition rate from state m to n 
and is generally dependent on the external force h. The 
time evolution operator 



■'van — ^nm 0~mn ^ ^ ^ink- 



(20) 



If the external force h(t) acting on the system is weak, 
Taylor expanding about h = 0, we get 

£mn(h) = (Co) mn + h(t)(C\) 

ra n • 

In this relation 
(A) ill n 



u nra , ~^nm 



k 



where a mn = [d n ^w mn ]h=o gives the relative change of 
rates. Note that the system is prepared in a NESS at 
h = characterized by the stationary distribution (p n )s- 
Then Eq. [6] can be expressed in the discrete notation as 

m,n 

= (A(t 2 )M(h)) (21) 
where the Agarwal term is 



(Pm)s n 

_ \ - (Pn)s 

/ / \ ^nm^nm 
„ (Pm)s 



E 



mn^mn • 



(22) 



Now we use the above relation to derive the velocity 
form of MFDR. We assume a displacement d mn associ- 
ated with each transition from state m to n. This has the 
property d mn = —d nm and gives a definition of velocity 
v(t) = J2 m ^ _ T m )d n -i lTn HI]. A generalized detailed 
balance in presence of the external force h 



(h) _ (0) 



exp[/3 hd ri 



(h) (0) 
leads to the following useful relation 

We also utilize the probability current 

J mn — P martin Pn^nm " Jn 



(23) 



(24) 



(25) 



Then from Eq. [22] we find the velocity form of Agwarwal 
term, 



Mm = 



{Pn)s 
(Pm): 

(Pn). 



~Wnrn Oi-nm 



E 



mn^mn 



^ , W nrn (a mn + f3dn m ) - ^2 ' 

n ^ Pm)s 

$ ^ 7 : ^nmdnm ^ 1 7 \ {J run) s^-mn 



where 



(Pm)s 
/3(u m - Z/ m ), 



E 



(26) 



(Pn)s 



nm^nm j 



(Pm)s 



(27) 



These relations lead to the velocity form of MFDR 
R A {t2-h) = ^ A m [e £o(t2 ~* l) p(ti)] rnll [/3(w I „ - ^ m )] 
= ^<A(t a )[t;(*i)-i/(*i)l>. (28) 

Note that Eq. [25] agrees with the results obtained in 
Rcf.s @, It- In particular, Ref. used the MFDR ex- 
pressed in terms of stochastic entropy of the system s 
(Eq. [5]) to obtain Eq. [35] They used the total stochastic 
entropy s to t = s + s m where s m is the stochastic entropy 
of the medium and showed 

dhS m (t) = }^5(t - T m )d m -X,m = v(t) = 0~n{t),mV m 
m m 

dhStot{t) = U(t) = y~] b~ n ( t ),mVm (29) 
in 

where v m and v m are given by Eq. 1271 

Note the equivalence of Eq. [35J with Eq.Q3] Indeed the 
analogy of v described here with the local mean veloc- 
ity \j(x,t)/p(x,t)] in the Langevin system becomes even 
more clear when we compare the steady state average 
(v s ) = J2 m (Pm.)s v m with (v s ) = J2 m (P™)s v ™ an d find 



,) =T^ J„ 



^ww = E(^< 

mn mn 



(v s ). (30) 



This relation is the same as Eq. [16] obtained for the 
Langevin system. 

For a velocity-response Eq.[25Jreadily leads us to Eq.[T5l 
already obtained in the context of Langevin dynamics. 
This completes one of the main achievements of this pa- 
per - the Agarwal formalism leads to the same form of 
velocity-MFDR for discrete master equation and contin- 
uum Langevin dynamics. 
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FIG. 1: Flashing ratchet model: The potential height 
switches between W = (off-state) and Wo (on-state). The 
asymmetry of the potential in the on-state is described by the 
inequality a ^ b. wi,2 denote transition rates between on and 
off states. 



VI. EINSTEIN RELATION 

Using the velocity MFDR (Eq. [Po|) we find the mobility 
in NESS 

poo poo 

H S = / drR v (r) = / dT(v( T )[v(0) - i/(0)]>. (31) 
Jo Jo 

On the other hand, the diffusion constant in an NESS 
having mean velocity (v s ) is 

/>oo 

D s = dT([v(r)-(v s )M0)-(v s )}). (32) 
Jo 

Thus the mobility fi s and diffusion constant D s at NESS 
do not satisfy the equilibrium Einstein relation, i.e., D s — 
Tji s = 1^0. The difference gives us the modification in 
the Einstein relation in terms of the violation integral 

poo 

I = D S -T^ S ^ dr [{v(t)v{0)) - (v s ) 2 } . (33) 
Jo 

Since the form of velocity-MFDR for Langevin equation 
(Eq. H5]) and master equation (Eq. |2"5)) are the same, 
we get the same modified Einstein relation for both the 
cases. 



VII. FLASHING RATCHET MODEL OF 
MOLECULAR MOTORS 

In this section we apply the concepts developed so 
far in this paper on a specific realization of the flash- 
ing ratchet model of molecular motors [38l - l40l l43l . |44| . 
In particular, we calculate the velocity-MFDR for this 
model and derive the violation integral of the correspond- 
ing non-equilibrium Einstein relation. 

A molecular motor, e.g., kinesin, moves along a poly- 
meric track, e.g., microtubule in a strongly fluctuating 
thermal environment utilizing intrinsic local assymmetry 
of the track and chemical energy provided by hydroly- 
sis of ATP to ADP and a phosphate. The binding and 
hydrolyzing of ATP changes the strength of interaction 




0.0002 0.0004 0.0006 0.0008 0.001 
r(s) 

FIG. 2: (color online) Velocity correlations and response: 
velocity response k B T R v (t), correlation function of the ve- 
locity with the local mean velocity (v(t)u(0)), and velocity 
auto correlation function (v(t)v(O)) as a function of time. 
The parameter-values used to obtain these curves are A = 
8nm, a/A = 0.1, D = 0.009 Jim's -1 , k B T = 4.2pN nm, 
Wo = 18.85 k B T and simulation time-step St = 1.8 x 10 -6 s. 
The transition rates are chosen to be equal with Wi = 
W2 = 3536 s" 1 . With these parameter values we find D s = 
0.0084 /im 2 s _1 , k B T fi s = 0.0057 ^m 2 s _1 , and the violation 
integral I = 0.0027 um s~ . The mean velocity in steady 
state is (v s ) — 2.04 um/s. 



of the motor with the polymeric track [431 ] . Thus a sim- 
ple two-state approximation of the dynamics of motor- 
proteins was proposed [43l . |44| where the motor encoun- 
ters a locally asymmetric but globally periodic potential, 
whose height switches between a large and a small value. 

We consider a flashing-ratchet model where the sys- 
tem swithces between two states, (1) on-state: stochas- 
tic motion in an asymmetric piece-wise linear potential, 
(2) off-state: simple one dimensional diffusion (Fig. [IJ. 
The probability distributions in the two states pi,2(£,0 
evolve by Q 

dtPi + dxji = u 2 p2 - ^ipi 

9tP2 + dx]2 = -U2P2 + LdlPl 

where ji = -Dd x pi and j 2 = -D{p 2 d x (W/T) + d x p 2 ], 
and denote the transition rates. In the on-state, 
the potential W(x) is periodic W(x) = W(x + A) with 
period A = (a + b). Within one period, W(x) = (Wo/a)x 
if < x < a, and W(x) = {W /b)(X - x) if a < x < A. 

We perform molecular dynamics simulations of a par- 
ticle moving under the influence of the above-mentioned 
ratchet potential in the presence of a Langevin heat bath. 
We use stochastic switching between the on and off states 
with a constant switching rate uj\ = uj2- From this sim- 
ulation, in Fig. [21 we plot the velocity-response func- 
tion kBTR v {r) and the related steady-state correlations 
(v(t)v(O)), (v(t)v(O)) (Eq. [T5J) . The parameter values 
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we use are enlisted in Fig. [2] and are typical of micro- 
tubule associated molecular motors J4f|. At long time, 
both (v(t)v(O)) and (w(r)^(O)} decorrelates to {v s } 2 . Wc 
utilize the correlation functions to determine the mobil- 
ity p Sl diffusion constant D s , and the violation integral 
/ = D s - k B Tfj, s (see Fig. EJ. 

We calculate the dependence of the steady-state mo- 
bility ksT fi s , diffusion constant D s and the violation 
integral / on the asymmetry parameter a = a/X (Fig. [3]) 
where a = 1/2 denotes the symmetric ratchet. This cal- 
culation leads us to the curious result that all the three 
quantities have minimum at a = 1/2. The steady-state 
diffusion constant D s in the flashing-ratchet is always 
suppressed (D s < D), and moves closer to the free dif- 
fusion D for the most asymmetric ratchet. Note that 
the violation integral quantifies the difference between 
NESS and equilibrium, with equilibrium requiring 1 = 0. 
The symmetric ratchet does not generate unidirectional 
motion, but the switching between the on and off states 
keeps the system out of equilibrium. Thus, though the 
violation integral reaches its minimum at a = 1/2 it re- 
mains 1^0. Setting switching rates loi,lo2 = would 
restore equilibrium with 1 = 0. The dependence of / on 
various models and parameter values at different NESS 
is yet to be fully understood. 



0.01 




FIG. 3: (color online) Flashing ratchet: diffusion constant 
D 3 , mobility fesT/i s and violation integral J as a function 
of asymmetry parameter a = a/\. All the other parameter 
values are same as in Fig [2] 



While calculation of all the other quantities from our 
simulations are straight-forward, u{t) demands a special 
mention. The local mean velocity v(x) is the stochastic 
particle velocity v averaged over the subset of trajectories 



passing through x. At steady state, this definition is 
the same as v s (x) = j s /Ps(z) where the mean current 
is constant everywhere: j s = p(v s ) with (v s ) the mean 
velocity at steady state and p = 1/A the mean density. In 
calculating (v(t)i/(0)), the local mean velocity at time t 
is obtained by identifying the value of v(x) corresponding 
to the position x visited by the particle at that instant. 



VIII. SUMMARY 

We have presented a unified derivation of modified fluc- 
tuation dissipation relations (MFDR) at non-equilibrium 
steady states (NESS) using the Agarwal formalism. Thus 
all the various versions of MFDR that we derived in 
this paper arc intrinsically equivalent to each other. Wc 
showed that the response function around any NESS can 
be expressed as a correlation between the observable and 
a variable conjugate to the external force with respect 
to the system's stochastic entropy production. For both 
a continuum Langevin and a discrete master equation 
system, we have shown that the non-equilibrium form of 
FDT involving velocity response can be expressed as an 
equilibrium one and an additive correction. The correc- 
tion in both these cases is a correlation function of the 
velocity with a local mean velocity. The resulting mod- 
ification of the Einstein's relation gives the violation in 
terms of a time integral over this additive correction. 

Using molecular dynamics simulations in presence 
of Langevin heat-bath, wc studied a flashing ratchet 
model within this framework and obtained the re- 
sponse function and velocity correlations in the steady 
state. We showed that the violation integral varies 
non-monotonically with the asymmetry parameter of the 
ratchet and reaches a non-zero minimum for the case 
of a symmetric ratchet. We plan to extend our study 
to other models of molecular motors stochastic 
particle-pumps pol |47| . polymer translocation dynam- 
ics [48| , and dynamics of self-propelled particles \4S^ . 
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